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Abstract. We report on a study of the anisotropy 
in the arrival direction of cosmic rays with a median 
energy per Cosmic Ray (CR) particle of ~ 14 TeV 
using data from the IceCube detector. IceCube is 
a neutrino observatory at the geographical South 
Pole, when completed it will comprise 80 strings 
plus 6 additional strings for the low energy array 
Deep Core. The strings are deployed in the deep 
ice between 1,450 and 2,450 meters depth, each 
string containing 60 optical sensors. The data used 
in this analysis were collected from April 2007 to 
March 2008 with 22 deployed strings. The data 
contain ~ 4.3 billion downward going muon events. 
A two-dimensional skymap is presented with an 
evidence of 0.06% large scale anisotropy. The energy 
dependence of this anisotropy at median energies per 
CR particle of 12 TeV and 126 TeV is also presented 
in this work. This anisotropy could arise from a 
number of possible effects; it could further enhance 
the understanding of the structure of the galactic 
magnetic field and possible cosmic ray sources. 
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I. Introduction 

The intensity of Galactic Cosmic Rays (GCRs) have 
been observed to show sidereal anisotropic variation on 
the order of 10~ 4 at energies in the range of 1-100 
TeV (HI, ||2l and 0). This anisotropy could arise from 
number of different combination of causes. One possible 
cause could be the Compton Getting (CG) effect. This 
effect was proposed in 1935 |4| predicting that CR 
anisotropy could arise from the movement of the solar 
system around the galactic center with the velocity of 
~ 220 kms^ 1 such that an excess of CR would be 
present in the direction of motion of the solar system 
while a deficit would appear in the opposite direction. 
Another possible effect (proposed by Nagashima et. 
al. 0) causing the excess in the anisotropy, which was 
referred to as "tail-in", originates from close to the tail 
of the heliosphere. While the deficit in the anisotropy, 
which was referred to as "loss-cone", originates from a 
magnetic cone shaped structure of the galactic field in 
the vicinity region. 

In this paper we present results on the observation of 
large scale cosmic ray anisotropy by IceCube. Previous 
experiments have published a 2-dimensional skymap of 
the northern hemisphere sky (0- (5)). This measure- 
ment presents the first 2-dimensional skymap for the 
southern hemisphere sky. In addition, we present the 



energy dependence of this anisotropy at median energies 
per CR particle of 12 TeV and 126 TeV. 

The outline of the paper is: the second section will 
describe the data used in this analysis, the analysis 
method and the challenges. The third section will discuss 
the results and the stability checks applied to the data. 
The fourth section will discuss the anisotropy energy 
dependence and the last section is the conclusion. 

II. Data Analysis 

The data used in this analysis are the downward going 
muons collected by the IceCube neutrino observatory 
comprising 22 strings. The data were collected from 
June 2007 to March 2008. The events used in this 
analysis are those reconstructed by an online one iter- 
ation Likelihood (LLH) based reconstruction algorithm. 
The events selected online require at least ten triggered 
optical sensors on at least three strings. The average rate 
of these events is ~ 240 Hz (approximately 40 % of the 
events at triggering level). Further selection criteria are 
applied to the data to ensure good quality and stable 
runs. The final data set consists of 4.3 x 10 9 events 
with a median angular resolution (angle between the 
reconstructed muon and the primary particle) of 3° and 
a median energy per CR particle of 14 TeV as simulated 
according to CORSIKA using SIBYLL hadronic 
interactions model and Horandel [8 | primary cosmic ray 
spectrum. 

In this analysis we are searching for a high precision 
anisotropy. The sidereal variation of the CR intensity 
is induced by the anisotropy in their arrival direction. 
However, it can also be caused by the detector exposure 
asymmetries, non-uniform time coverage, diurnal and 
seasonal variation of the atmospheric temperature. Apart 
from these effects the remaining variations can only be 
of galactic origin. 

Due to the unique location of IceCube at the South 
Pole the detector observes the sky uniformly. This is not 
the case for all the previous experiments searching for 
large scale cosmic ray anisotropy (e.g. 0, 0, 0). 
Due to their locations they need a whole solar day to 
scan the entire sky. As a result they need to eliminate the 
diurnal and seasonal variations using various techniques. 
For IceCube the sidereal variation is not affected by 
diurnal variations because the whole sky is fully visible 
to the detector at any given time and because there 
is only one day and one night per year. In addition, 
although the seasonal variation of muon event rate is 
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on the order of 20% the variation is slow and does not 
affect the daily muon intensity significantly. 

The remaining challenge for this analysis is account- 
ing for the detector asymmetry, and unequal time cover- 
age in the data due to the detector run selection. To illus- 
trate the detector asymmetry Figure Q] shows the IceCube 
22 string geometrical configuration. This geometrical 
asymmetry results in a preferred reconstructed muon 
direction since the muons would pass by more strings 
in one direction in the detector compared to another. 
The combination of detector event asymmetry with a 
non-uniform time coverage would induce an azimuthal 
asymmetry and consequently artificial anisotropy of the 
arrival direction of cosmic rays. This asymmetry is 
corrected for by normalizing the azimuthal distribution. 

Figure |2] shows the azimuthal distribution for the 
whole data set. It displays the number of events vs. 
the azimuth of the arrival direction of the primary 
CR particle. Note that the asymmetry in the azimuthal 
distribution due to detector geometry is modeled well by 
simulation. 

To correct for the detector azimuthal asymmetry we 
apply an azimuthal normalization. The azimuthal distri- 
bution is parameterized by N, nu and n, where N is the 
total number of bins, is the number of events per bin 
and n is the average number of events, n = -k Y^iLi n i- 
n is denoted by the horizontal red line in Figure [2] The 
azimuthal normalization is applied by weighting each 
event by — for that event. 

In addition to the azimuthal asymmetry we also 
observe a non uniform zenith angle distribution (more 
events arrive from the zenith than from the horizon). 
Due to this declination dependence, the sky is divided 
into four declination bands such that the data is approxi- 
mately equally distributed among the declination bands. 
For each band the azimuthal distribution is normalized 
for the whole year. The relative intensity for each bin in 
the 2-dimensional skymap is then calculated by dividing 
the number of events per bin by the total number of 
events per declination. 

III. Results: 

Figure [5] shows the southern hemisphere skymap 
for well reconstructed downward going muons for the 
IceCube 22-Strings data set. The skymap is plotted in 
equatorial coordinates. The color scale represents the 
relative intensity of the rate for each bin per declination 
band where each bin rate is calculated by dividing the 
number of events for that bin over the average number 
of events for that bin's declination band. The plot shows 
a large scale anisotropy in the arrival direction of cosmic 
rays. The amplitude and the phase of this anisotropy is 
determined by projecting the 2-dimensional skymap in 
Right Ascension (RA) as shown in Figure [4] Figure [4] 
shows the relative intensity vs. the RA. The data are 
shown in points with their error bars. The fit is the 
first and second-order harmonic function in the form 
of Yn=i{ A i x cos (* (( RA ) - <f>i))) + B where A i is 
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Fig. 1: The IceCube detector configuration. The filled 
green circles are the positions of IceCube strings and 
the filled blue circles display the position of the IceTop 
tanks. 
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Fig. 2: The azimuthal distribution for the whole data set. 
This plot shows the number of events vs. the azimuth 
of the arrival direction of the primary CR particle. The 
horizontal red line in the average number of events for 
the distribution. 



the amplitude and fa is the phase and B is a constant. 
The fit Ai, fa and \ 2 /ndof for the first and second 
harmonic fit are listed in table H] The significance of 
the 2-dimensional skymap is shown in Figure [5] The 
significance is calculated for each bin from the average 
number of events for that bin's declination band. Note 
that the significance of several bins in the excess region 
is greater than 4<r and in some bins in the deficit region 
is smaller the — 4er. 

To check for the stability of the measured large scale 
anisotropy we performed a number of checks with the 
data set. One check was applied by dividing the data into 
two sets where one set contains sub-runs with an even 
index number and the other set contains sub-runs with 
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Fig. 3: The IceCube skymap in equatorial coordinates 
(Declination (Dec) vs. Right Ascension (RA)). The color 
scale is the relative intensity. 



In addition, stability checks are applied to test for 
daily and seasonal variation effects. To test for the daily 
variation effect the data were divided in two sets: The 
first set contains sub-runs with rate values greater than 
the mean rate value for that sub-run's day. The second 
set contains sub-runs with rate values less than the 
mean rate value for that sub-run's day. Furthermore, 
to test for the seasonal variation effect the data set is 
divided in two sets: The first set holds the winter month's 
sub-runs (June-Oct.). The second set holds the summer 
month's sub-runs (Nov.-March). For both tests we see 
no significant changes in the value of the anisotropy of 
the two data sets. 
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Fig. 4: The 1-dimensional projection of the IceCube 2- 
dimensional skymap. The line is the first and second 
harmonic function fit. 
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TABLE I: The first and second harmonic fit amplitude, 
phase, and x 2 /ndof. 
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Fig. 5: The IceCube significance skymap in equatorial 
coordinates (Dec vs. RA). The color scale is the signif- 
icance. 



an odd index number (a sub-run on average contains 
events collected for ~ 20 minutes at a time). Another 
check is applied by dividing the data into two sets each 
set contains half the sub-runs selected randomly. The 
results for both tests were consistent. 



IV. Energy Dependence: 

In order to better understand the possible nature of 
the anisotropy we have searched for energy dependent 
effects using our data. To determine the energy depen- 
dence for the signal we divided the data into two energy 
bins. To accomplish that and to ensure constant energy 
distribution along our sky both the number of sensors 
triggered by the event and the zenith angle of the event 
are used for the energy bands selection. The first energy 
bin contains 3.8 x 10 9 events with a median per CR 
particle of 12.6 TeV and 90% of the events between 2 
and 158 TeV. The second energy bin contains 9.6 x 10 8 
events with a median energy per CR particle of 126 TeV 
and 90% of the events between 10 TeV and 1 PeV. Each 
2-dimensional skymap is projected to 1-dimensional 
variations in RA. In comparison to previous experiments 
the 1-dimensional RA distribution is fitted to a first 
harmonic fit. The first harmonic fit amplitude and phase 
for the first energy band are A x = (7.3 ± 0.3) x 10~ 4 
and (pi — 63.4° ± 2.6°, while the amplitude and phase 
for the first energy band are A\ = (2.9 ± 0.6) x 10~ 4 
and (j>i — 93.2° ± 12°. Figure [6] shows the amplitude vs. 
energy of this analysis (in filled circles) in comparison 
to previous experiments (In empty squares). Note that 
the amplitude in this analysis shows a decrease of the 
harmonic amplitude value at the higher energies for the 
energy ranges of 10-100 TeV. 

V. Conclusion 

In this analysis we present the first 2-dimensional 
skymap for the southern hemisphere of 4.3 billion cos- 
mic rays with a median angular resolution of 3° and a 
median energy per CR particle of 14 TeV as observed 
by IceCube. A large cosmic ray anisotropy with a first 
harmonic vector amplitude of A\ — (6.4 ± 0.2) x 10 -4 
and a phase of (f>i = 66.4° ± 2.6° is observed. The 
significance of some bins in the excess and the deficit 
regions were found to be > |4er|. This anisotropy is an 
extension of previously measured large scale anisotropy 
at the northern hemisphere reported by multiple experi- 
ments (0- 0). 

In addition, we report on the anisotropy energy de- 
pendence. We report the amplitude of the first harmonic 
vector of the anisotropy for the two energy bands. 
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Fig. 6: The amplitude vs. energy in eV. The filled circles 
markers are the result of this analysis and the empty 
square markers are the result from previous experiments 

(0, AH, Col, CD, G2, US, d, G3, EE 03, ED, 
El, Eol, ED, E2, ED, El, E3). 



The first energy band with a median energy per CR 
particle of 12.6 TeV the amplitude is found to be 
Ai = (7.3±0.3) x 1CP 4 . The second energy band with a 
median energy per CR particle of 126 TeV the amplitude 
is found to be A 1 = (2.9 ± 0.6) x 1CT 4 . The amplitude 
energy dependence is found to follow a decreasing trend 
with energy. 
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